Particle flows in a dc discharge in laboratory and microgravity conditions 
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We describe a series of experiments on dust particles flows in a positive column of a horizontal 
dc discharge operating in laboratory and microgravity conditions. The main observation is that 
the particle flow velocities in laboratory experiments are systematically higher than in microgravity 
experiments, for otherwise identical discharge conditions. The paper provides an explanation for 
this interesting and unexpected observation. The explanation is based on a physical model, which 
properly takes into account main plasma-particle interaction mechanisms relevant to the described 
experimental study. Comparison of experimentally measured particle velocities and those calculated 
using the proposed model demonstrates reasonable agreement, both in laboratory and microgravity 
conditions, in the entire range of discharge parameters investigated. 



I. INTRODUCTION 

Complex (dusty) plasmas constitute an interdisci- 
plinary research field with relations to space and as- 
trophysical topics |l|, industrial plasma applications [3, 
plasma fusion oriented research [3|, physics of strongly 
coupling phenomena [1], and soft condensed matter 0,[al- 
The problem of plasma-particle interactions is central to 
this field, because these interactions govern practically all 
phenomena that can be observed and investigated [3-Q . 
Of particular importance are interactions that affect par- 
ticle charging and screening [9|, |lO[ , momentum transfer 
between different complex plasma components [111, ex- 
ternal and internal forces acting on the particles [9|. 

The focus of this paper is on a situation, where dif- 
ferent complex plasma components are allowed to drift 
relative to each other, in an electric field of a positive col- 
umn of a direct current (dc) discharge. Particles are in- 
jected into a horizontally mounted dc discharge tube and 
are transported along the tube by various forces, most 
important of which are the electrical, ion drag, and neu- 
tral drag forces. The particle velocities can be relatively 
easily measured in a wide range of discharge parameters 
and they provide important information on various basic 
complex plasma properties, in particular particle charges 
and forces they are acted upon. A series of experiments 
that we describe here is performed using the same exper- 
imental setup in laboratory and under microgravity con- 
ditions (during the parabolic flights). The main obser- 
vation is that particle velocities are systematically higher 
in laboratory as compared to microgravity conditions at 
otherwise identical discharge conditions. This finding is 
naturally a puzzling one, since the plasma parameters 
(and hence plasma-particle interaction mechanisms) are 
not expected to depend on the presence or absence of 
gravity. The main purpose of the present paper is to 
provide a convincing explanation for this unexpected ob- 
servation. 

We analyze in detail the specifics of particle fiows in a 
positive column of a horizontal dc discharge in ground- 
based (laboratory) and microgravity conditions. In doing 



so we put forward a physical model, which is believed 
to correctly describe main plasma-particle interactions 
relevant to the present experiment. The issue of parti- 
cle charging in a flowing collisional plasma is the main 
constituent of the model. We show that in laboratory, 
the force of gravity shifts particles downwards from the 
vicinity of the tube axis. As a result, parameters of the 
plasmas surrounding the particles in laboratory and mi- 
crogravity conditions are somewhat different. In partic- 
ular, particles in laboratory experiments are located in a 
region with higher ion drift velocity. Higher ion drift ve- 
locities result in higher particle charges and higher elec- 
trical forces acting on the particles. This is the main 
qualitative explanation behind the difference in particle 
velocities measured in laboratory and microgravity con- 
ditions. We finally demonstrate reasonable quantitative 
agreement between experimental measurements and the 
analytical theory proposed in this paper. 

The paper is organized as follows. In Section |ll] we 
describe the experimental apparatus, experimental pro- 
cedures, and main observations. Section IIIII introduces 
the model of particle charging relevant to the present 
experimental conditions. In Section IIVI we summarize 
main forces acting on the particles and in Section FVl we 
consider the force balance that determines the particle 
motion, separately for laboratory and microgravity con- 
ditions. Results from this analysis are discussed in Sec- 
tion IVII This is followed by conclusion in Section IVIII 



II. EXPERIMENT 

The experiments are carried out in the Plasma Kristall- 
4 (PK-4) facility, scheduled for operation onboard the 
International Space Station (ISS) after 2014 [Ulil. It 
represents a dc discharge generated in a U-shaped glass 
tube of 3 cm in diameter. The main part of the tube 
containing the positive column of the discharge is 35 cm 
long. In the described experiments the horizontal con- 
figuration of the tube is employed (see Fig. [T]). Complex 
plasma is formed by injecting /im-size grains into the 
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FIG. 1: Sketch of the PK-4 facility in the horizontal configu- 
ration. 



discharge. The particles can be observed and recorded 
by video cameras at a required position inside the tube. 
Normally, individual particle resolution is possible. The 
sketch of the PK-4 experimental setup is shown in Fig. [T] 
For a cornprehensive description of this project see e.g. 
Refs. [ll,|ll. 

The PK-4 facility is particularly suited to study dif- 
ferent kinds of flow-related phenomena in complex plas- 
mas. Previous experiments in laboratory and micro- 
gravity (parabolic flights) conditions have already pro- 
vided a wealth of useful information, in particular, re- 
garding particle charging [lJ-[3, the ion drag force [l3] , 
particle-particle interactions jlSl [lOf, linear and non- 
linear wave s |20|. 2l| , structural properties of the particle 
component [2^[23|, and other important phenomena [2J]. 
Experiments reported here correspond to a series of ba- 
sic studies, aiming at checking the reliability of the PK-4 
setup under microgravity conditions 

The experiments have been performed in laboratory 
conditions and under microgravity, during the parabolic 
flights campaign in the fall 2012. For these experi- 
ments neon gas in the pressure range between ^ 30 
and ^ 90 Pa is used, the discharge current is fixed at 
1 mA. Spherical melamine-formaldehyde particles (mass 
density is p ~ 1.51 g/cm^) are injected into the discharge 
tube. Being injected into the plasma, the particles be- 
come charged negatively and are transported through the 
positive column of the discharge by the combination of 
various forces. The dominant forces in the horizontal 
direction are the electrical, ion drag and neutral drag 
forces p^| - [la |. The particles are illuminated by a laser 
sheet and their motion is recorded by a video-camera, 
situated near the center of the tube. The video-camera 
has a field of view of 22.1 x 16.6 mm^ and operates at 
a frame rate 35 fps. Figure [5] shows typical snapshots of 
drifting particle clouds in microgravity (left panel) and 
laboratory conditions (right panel) . 

The velocities of the particles are estimated by measur- 
ing their track lengths in single images (snapshots). Each 
experimental point corresponds to the average value of 
track lengths selected randomly in different images for 
a particular set of experimental parameters (neutral gas 



pressure, particle size, particle dispenser). The procedure 
is essentially the same as used previously [ij, [l^ • 

The main observations are as follows: (i) For any par- 
ticular particle size the particle drift velocities decrease 
monotonically with increasing pressure; (ii) For a given 
particle size and pressure, the drift velocities in labora- 
tory experiments are systematically higher than in mi- 
crogravity conditions; (ii) The difference between drift 
velocities in laboratory and microgravity conditions in- 
creases with increasing the particle size and with lowering 
neutral gas pressure. 

To proceed with quantitative explanation of these ob- 
servations, the knowledge of main plasma parameters is 
required. The plasma density, electron temperature, and 
axial electric field were measured previously in the ab- 
sence of particles [29| using a single Langmuir probe of 
~ 4 mm length and ~ 30 /zm in diameter, insulated by 
a glass holder [SO]. The linear fits to the experimen- 
tal results applicable in the pressure range investigated 
were given in Ref. [l5|. Since a relatively small num- 
ber of particles is injected in each experimental run and 
the particle cloud extent in the radial direction is con- 
siderably smaller than the tube radius (see Fig. [2]), we 
neglect all kinds of collective effects as well as the effect 
of the particles on plasma parameters. The analysis is 
thus performed in the individual particle approximation. 



III. MODEL OF PARTICLE CHARGING 

The purpose of this Section is to present a simple 
heuristic model of particle charging in collisional plas- 
mas with flowing ions, which can be used to explain the 
observations summarized above. Two important effects 
should be accounted for in such a model: Ion-neutral 
collisions and relative drift between the ions and the par- 
ticle. Although, the effect of ion-neutral collisions on the 
particle charging has received considerable attention |3ll - 
|37| (for a review see Ref. [l^l ) , most of the studies have 
concentrated on isotropic conditions. To the best of our 
knowledge there is very few papers, which reports re- 
sults for ion collection by a sphere in a drifting collisional 
plasma [3^, |33. Numerical results from Ref. [33 are 
more comprehensive from the point of view of compari- 
son, and thus they will be used to check the reliability of 
the approximation we put forward here. 

The reduced particle surface potential z — e\(j)s\/Te 
(and thus the particle charge Q ~ ^ga) is set by the bal- 
ance between the electron and ion ffuxes absorbed on the 
particle surface (here 0s is the particle surface potential, 
a is the particle radius, T^ is the electron temperature, 
and e is the elementary charge). In order to simplify the 
comparison with Ref. [39| we use the following normaliza- 
tion for the fluxes: The fluxes are in u nits of iTrnja^Cs, 
where rii is the ion density and Cg = ^T^jmi is the ion 
sound speed. In the pressure range investigated, electron- 
neutral collisions are not important and the conventional 




FIG. 2: Example of the particle flow in microgravity (left panel) and laboratory (right panel) experiments. Both snapshots 
correspond to the particle size 2a — 3.43 /im and neutral gas pressure p — 50 Pa. Particles are drifting to the left. The particle 
flow velocity in laboratopry conditions is almost a factor of two higher (note longer track lengths). Note also that in laboratory 
conditions the particle cloud is considerably shifted (downwards) from the tube axis, which corresponds to the centers of the 
images. Waves visible in the lower part of the right panel are related to an increase in the total electric fleld (radial field 
is required to levitate the particles in laborator y c onditions). These waves are excited by sufficiently fast ion fiows and the 
corresponding ion-dust streaming instability [25l - [2a |. The field of view is 22.1 x 16.6 mm (1600 x 1200 pixels ). The exposure 
time is 28 ms. 



OML expression for the electron flux [40| is applicable, 



je = V?7ii/27rmeexp (-z) . 



(1) 



In order to derive an approximation for the ion flux 
we consider separately two regimes. In the first - weakly 
collisional (WC) regime - a reasonable approximation for 
the ion flux is a linear superposition of the coUisionless 
(OML) and collisional contributions [3l,|3l|. The OML 
expression for the shifted Maxwellian distribution reads 
in the present notation as 
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where £_ = M-^t/2 is the measure of the ion drift velocity 
expressed in terms of the Mach number, M — u/Cs, and 
T = Te/Ti is the electron-to-ion temperature ratio. The 
collisional contribution to the ion flux in the isotropic sit- 
uation can be roughly estimated (in dimensional units) 
as A Jf°'' ~ (47r/3)i?on,;i', where Rq denotes the radius of 
a sphere inside which the ion-particle interaction is suf- 
ficiently strong [32, |41| and v is the effective ion-neutral 
collision frequency. For the Debye-Hiickel (Yukawa) po- 
tential a simple formula of the type Rq ~ X\n{l + Rc/X), 
where A is the plasma screening length and Rq = \Q\e/Ti 
is the (ion) Coulomb radius, would describe adequately 
the respective limits of weak and strong ion-particle cou- 
pling and provide a smooth transition between them [42] . 
The remaining step is to estimate how the presence of the 
ion flow modifies this collisional contribution. 

In the spirit of Ref. [ij] we define the effective Coulomb 
radius and plasma screening length in flowing plasma as 



follows. The Coulomb radius is Re 



\Q\e 



SO that 



it reduces to the standard definition in isotropic situ- 
ation (^ = 0) and decreases as the ion kinetic energy 
(associated with the direct ion flow) increases. The effec- 
tive plasma screening behaves as A ~ Aoe/i/l + fiO^, 
where /(^) = (1 + 2^^)~^ is an adjusting function and 
Ade = \/Te/4:Tre'^ne is the electron Debye radius. This 
ensures that A tends to Adj: = y/Ti/Am?ni in the 
isotropic regime, approaches Aoe for highly suprather- 
mal flows, and provides smooth transition between these 
limits [4^,|4J1- With the use of these approximations, the 
collisional contribution to the ion flux can be written in 
the present (dimensionless) notation as 
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where Uc = va/Cs is the reduced ion-neutral collision fre- 
quency (ion coUisionality index) and L,, = Xj^^/a. The 
total flux that the particle collects in the weakly colli- 
sional regime is simply jwc — Ji^^^ + Aj,™". 

In the highly collisional (HC) regime we make use of 
the problem of collisional ion flow around a Coulomb 
sphere, which allows for an analytical solution [45|. In 
the present notation the result can be written as: 
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Finally, the interpolation formula, which is applicable 
in the entire range of ion coUisionality is 



Ji — Uwc + .?HCJ 



-T.x-l/7 



(5) 
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FIG. 3: (Color) Normalized ion flux to the particle surface ji 
(top panel) and reduced particle charge z (bottom panel) vs. 
the reduced ion-neutral collision frequency ujc for several ion 
drift velocities. Symbols correspond to the numerical results 
from Ref. [3y]; solid curves are computed using the present 
approximation [Eq. [S]. At the left border of the top panel 
the symbols and curves from top to bottom correspond to 
M = 0.0 (red), M = 0.2 (blue), M = 0.4 (green), M = 0.8 
(magenta), M — 1.6 (cyan), and M = 3.2 (black). Similarly, 
at the left border of the bottom panel the curves and symbols 
from top to bottom are for M = 1.6, M = 0.8, M = 0.4, 
and M = 0.0, respectively. The color schemes are identical in 
the top and bottom panels (bottom panel shows less results 
for clarity). The plasma parameters used in simulations and 
analytical estimates are: H+ ions, Le = 20 and r — 10. 



approach developed, we performed calculations for the 
parameter regime studied by Haakonsen and Hutchin- 
son [331 ■ The results of this comparison are shown in 
Fig. El The top panel demonstrates that the proposed ap- 
proximation provides fairly good agreement with the nu- 
merical results in the regimes of weak and strong ion col- 
lisionality, but overestimates the ion fluxes in the vicinity 
of their peaks (except the case of highly supersonic ion 
drift, where the ion flux is underestimated). The bot- 
tom panel shows the corresponding values of the reduced 
particle charge, calculated from z = — ^ \n(27Tmejf /Trii). 
The agreement between the numerical simulations and 
the present analytical approximation is rather good in 
the entire regime of ion coUisionality. We should men- 
tion that the numerical results from Ref. [3y| are for a 
particular scenario of ion flow generation. Namely, in 
simulations ion drift is generated by ion-neutral charge 
exchange collisions in the presence of drifting neutral 
background. We expect, however, that the approxima- 
tion described above is not very sensitive to such details 
and the essential physics is captured. 

The important result, documented in Fig. [31 is that 
in the regime of weak and moderate ion coUisionality 
{ujc ^ 1) the ion flux collected by the particle decreases 
with the ion flow velocity. This decrease is especially 
pronounced at moderate coUisionalities, where significant 
coUisional enhancement of the ion flux is observed in the 
isotropic case. This enhancement is weakened as the 
ion fiow velocity increases. The physical reason is clear: 
Higher flow velocities imply higher ion kinetic energies, 
which makes the region of strong ion-particle interaction 
narrower, and, as a result, reduces the collisional contri- 
bution to the ion flux. Regarding the particle charge, its 
absolute magnitude increases with the ion flow velocity 
in the considered regime of weak and moderate ion coUi- 
sionality (this is of course the consequence of decreasing 
the flux of positive ions to the particle surface). Note, 
that the trend will reverse in the regime of strong ion 
coUisionality, which is, however, of no interest for the 
present study. 

The ion-neutral collision frequency v has been treated 
as an independent parameter in this Section. In fact, 
however, it is generally dependent on the ion drift ve- 
locity. To simplify the calculations we will employ the 
assumption of the constant ion- neutral collision cross sec- 
tion below. The collisional frequency is approximated as 



% 



(6) 



where 7 is generally an adjustable parameter [33|. The 
simplest choice 7 = 1 has been shown to yield reasonable 
results in many cases (Sg and we adopt it here, although 
it should be realized that treating 7 as a free parameter 
can naturally provide better accuracy [la, [3j] . 

The reduced surface potential of the particle is then ob- 
tained by equating the electron and ion fluxes [Eqs. ([T]) 
and ©J respectively]. As a check of the reliability of the 



so that 1/ ~ n„avTi ^^ ^^^ isotropic regime and ;/ ~ n„crw 
in the highly anisotropic regime. The cross section a ~ 
1 X 10""'^"' is chosen, which would reasonably describe the 
mobility of Ne"*" ions in neon in the near-thermal drift 
regime (see Appendix). 

The last important point to mention is that ion-neutral 
collisions are not the only mechanism leading to the en- 
hanced ion collection. Essentially the same effect can 
be associated with electron impact ionization events in 



the vicinity of the particle I4l|. Both mechanisms can 
be added in a simple superposition and their cumulative 
effect is Aj™"^ = (1 + /C)j™", where /C = vi/v is just 
the ratio of ionization and ion-neutral collision frequen- 
cies. This ratio is a sharply increasing function of the 
electron temperature and exceeds unity for Tg <^ 6.5 eV 
in neon plasmas [41| . Since the electron temperatures up 
to Te ~ 8 eV have been measured in the PK-4 facility in 
the regime of low neutral gas pressures [l^ , the effect of 
ionization should not be neglected in the present study. 
We, therefore, take it into account, when evaluating the 
particle charge. 



IV. MAIN FORCES ON THE PARTICLES 

The main forces affecting charged grains in complex 
plasmas subject to an external dc electric field are the 
electrical force and the forces associated with the mo- 
mentum transfer from drifting ions and electrons - ion 
and electron drag forces. The electrical force is 



Fci — QE, 



(7) 



where the particle charge Q implicitly depends on the 
electric field via the ion drift velocity u. The ion drift ve- 
locity and the electric field arc related via the ion mobility 
coefficient, u = /iE. In the Frost approximation [4g| we 
have 



u ~ a{E/p) [I -f- b{E/p)] 



'1/2 



(8) 



where the coefficients are a ~ 1.4 x 10^ dyne/statvolt/s 
and 6 ~ 1.6 x 10^ dyne/statvolt/cm (in CGS units). The 
ion flow, caused by the electric field, naturally results in 
the momentum transfer from the ions drifting relative to 
the particle component. The resulting force - the ion 
drag force [47H49j - can be written as 
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(9) 



where Uid is the characteristic momentum transfer fre- 
quency in ion-particle collisions. The relative velocity is 
very close to the actual ion drift velocity since the parti- 
cles are very massive and their motion can be neglected 
in comparison with that of the ions. For the momentum 
transfer frequency we use the result derived in Ref. [43l | 
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where In A is the Coulomb logarithm. This expression 
is applicable for arbitrary ion drift velocity, but does not 
take into account the effect of ion- neutral collisions on the 



ion drag force [38|, |5CH52| . The latter is relatively small as 
long as the ion mean free path does not become very short 
compared to the plasma screening length [50| . Therefore, 
for the present purposes we can neg lect coUisional effects. 
The Coulomb logarithm reads |43j 



lnA~ln 



Rc + X 



Rc+a 



(11) 



where the Coulomb radius Rq and the effective plasma 
screening length A depend on the ion drift velocity as 
discussed in Sec. IIIII 

Electric fields can also generate electron drifts (in par- 
ticular, this occurs in the axial electric field in the positive 
column of dc discharges, which is relevant to the present 
study). However, the corresponding electron drag force 
is usually much smaller than the ion drag force, provided 
the electron-to-ion temperature ratio is high j53| . Thus, 
the electron drag force is neglected in thepresent study. 

The neutral drag force (Epstein drag [5J]) acts in the 
direction opposite to the particle motion. It can be writ- 
ten as 



-VndmnVd, 



(12) 



where r?T,„ is the neutral mass, Vd is the particle ve- 
locity relative to the stationary background of neu- 
trals (gas flows are absent in these experiments), and 
Vnd = (8v27r/3)(5a^n„WT„ is the momentum transfer fre- 
quency in particle-neutral collisions. The numerical fac- 
tor (5 = 1 + ^~1.4 corresponding to diffuse scatter- 
ing with complete accommodation of neutrals from the 
particle surface is used, in agreement with the reported 
experimental results on /im-size melamine formaldehyde 
spheres [SSj . 

Finally, in laboratory experiments particles are natu- 
rally affected by the force of gravity, Fg = rridg, where 
g ~ 980 cm/s^ is the gravitational acceleration. 



V. FORCE BALANCE IN LABORATORY AND 
MICROGRAVITY CONDITIONS 

The qualitative difference between the balance of forces 
acting on the particles in laboratory and migrogravity 
conditions is illustrated in Fig. |31 Under microgravity 
the particle cloud is located centrally (near the tube axis) 
and the particles feel only the horizontal electric field 
E\ I , which also determines the ion drift velocity. In con- 
trast, in laboratory conditions the particle cloud is some- 
what shifted downwards from the center (see Fig. ^. 
This is obviously because some radial electric field is 
required, in order for the corresponding electrical force 
pointing upwards can balance for the particle gravity. 
At this position the particles feel the total electric field, 
which is a sum of the horizontal and radial components, 

E-- 
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The total ion drift velocity at the posi- 
tion of particles is now determined by E, not by E\\ as in 
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FIG. 4: (Color online) Sketch of the forces acting on a single particle in microgravity (a) and laboratory (b) conditions. The 
longitudinal discharge electric field is pointing to the right. Negatively charged particle moves to the left. 



microgravity. Since E > E\\ the ion velocity is effectively 
higher for particles in ground-based conditions than un- 
der microgravity. As has been discussed in Section IIIIl 
higher drift velocities imply higher particle charges and 
therefore higher values of the electrical force. This merely 
explains why the particle drift velocity should be higher 
in ground based conditions than in microgravity for the 
same discharge conditions. (Note that the ratio of the ion 
drag force to the electrical force remains almost constant 
for subthermal ion drifts and then decreases rapidly with 
ion velocity in the suprathermal regime |43l |. and this 
can be an additional factor contributing to the differ- 
ence in particle velocities) . The discussed effect is clearly 
more pronounced for bigger particles, since bigger parti- 
cles imply higher E±^ and hence larger difference between 
E and E\\. It also should disappear at sufficiently high 
pressures, such that ion drifts are subthermal and the 
particle charge is independent of the ion drift velocity. 

Under microgravity conditions, the particles drift ve- 
locity can be simply calculated from 
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(13) 



Probe measurements have evidenced that the axial elec- 
tric field is practically independent of pressure and we 
therefore take E\\ ~ 2.1 V/cm as fixed. Then the calcu- 
lations proceed as follows: (i) The ion drift velocity, as a 
function of pressure, is calculated from Eq. ([5]); (ii) The 
ion- neutral collision frequency is calculated from Eq. ([B]) ; 
(iii) The particle charge, as a function of pressure, is eval- 
uated from the model described in Sec. IIIIl (talking into 
account the dependence of relevant parameters on pres- 
sure, ion drift velocity, and collisional frequency); (iv) 
The ion drag force is calculated from Eqs. (|9l)- (fTT|) : (v) 
Solution of Eq. [13] then yields the dependence of particle 
drift velocity on pressure. 

Regarding laboratory experiments, the calculations are 
more involved in this case, because the calculation of the 
particle charge and the ion drag force should now be cou- 
pled self- consistently to the force balance condition in the 
vertical direction 
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In practice. Equation (|14p along with the model equa- 
tions for particle charge and ion drag force are solved to 
determine E±^ and the total electric field E. This fixes 
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FIG. 5: (Color online) Particle drift velocity as a function 
of the neutral gas pressure for particles with diameter 2a = 
1.2 /im. Symbols correspond to the previous (laboratory) 
experimental results [ij, [I3 , solid curve is obtained using the 
model of the present paper. 



the total ion drift velocity, collisional frequency, particle 
charge, and the ion drag force so that the particle veloc- 
ity can be evaluated from the force balance condition in 
the horizontal direction [Eq. ([T3|)]. 

Results from these calculations will be presented in the 
next Section. 



VI. RESULTS 

We start with presenting results for the case of small 
particles of diameter 2a ~ 1.2 ^m. For such small parti- 
cles gravity plays almost no role since the radial electric 
field E± required for the force balance is several times 
smaller than the axial electric field E\\ [l 



II 



For this 



reason, experiments under microgravity conditions have 
not been performed for this particle size. Figure [S] shows 
the comparison between experimental par ticle velocities 
measured previously in laboratory [ij, [l5| and theoreti- 
cal calculation using the model presented above (only one 
curve is shown since the curves for laboratory and micro- 
gravity conditions are essentially coinciding) . The agree- 
ment is convincing, indicating that the necessary ingre- 
dients have been properly incorporated into the model. 

Next, we consider larger particles of 2a ~ 2.55 /xm in 
diameter. The results of experiments and calculations 
arc presented in Fig. [5} The solid (open) circles corre- 
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FIG. 6: (Color online) Particle drift velocity as a function of 
the neutral gas pressure for particles with diameter 2a = 2.55 
fim. Open (red) circles correspond to the previous laboratory 
experimental results [Ij, ]1M- Solid (red) circles correspond 
to the present laboratory experiment. Solid (blue) squares 
are the measurements made under microgravity conditions. 
Upper (red) and lower (blue) curves are obtained using the 
model of the present paper for laboratory and microgravity 
conditions, respectively. 



spond to the present (previous [l4l [l5| ) experiments in 
laboratory. Solid squares are obtained from the analysis 
of parabolic flights experiments. The upper (red) curve is 
calculated taking into account gravity. The lower (blue) 
curve corresponds to calculations for microgravity con- 
ditions. Analytical results are in reasonable agreement 
with the experimental measurements and correctly re- 
produce the tendency for particles to drift faster in lab- 
oratory conditions. 

Finally, we consider the largest particle size investi- 
gated, 2a = 3.43 /xm. For these particles two different 
dispensers have been used to produce particle injection 
into the discharge. One is conventional shake dispenser 
(SD) and the other is the gas-jet dispenser (GJD). It is 
observed that particles drift faster when GJD is used, 
the relative difference can amount to ^ 50% ^3- Nev- 
ertheless, clear difference between particle flow velocities 
in laboratory and microgravity is still observed. This is 
illustrated in Fig. jT] Here the solid circles correspond 
to the laboratory experiments, while solid squares are 
from microgravity experiments. In each case, upper sym- 
bol for a given pressure corresponds to the use of GJD. 
The upper (red) curve is calculated taking into account 
gravity. The lower (blue) curve corresponds to calcu- 
lations for microgravity conditions. We see again that 
analytical results are in reasonable agreement with the 
experimental measurements. Comparing Figures IHl and 
[7] we also see that the difference in particle drift veloc- 
ities increases with the particle size and with lowering 
the neutral gas pressure. Thus the proposed theoretical 
model adequately describes the main experimental ob- 
servations. 

To conclude this section we comment on the sources of 
scattering of experimental points in Figs. |5]~|71 The first 
source is related to uncertainties in measured velocities, 
which are normally estimated as 10—15% for the pressure 



FIG. 7: (Color online) Particle drift velocity as a function of 
the neutral gas pressure for particles with diameter 2a = 3.43 
/im. Solid (red) circles correspond to the present laboratory 
experiment. Solid (blue) squares are the measurements made 
under microgravity conditions. Upper (red) curve and lower 
(blue) curves are obtained using the model of the present pa- 
per for laboratory and microgravity conditions, respectively. 



range investigated [T^ [T^ . The second possible source 
is related to the effect of particles on discharge parame- 
ters. The number of injected particles can vary from one 
experiment to another. Larger numbers of injected par- 
ticles imply larger modifications of surrounding plasma 
parameters. Although in the present experiment care has 
been taken to avoid high particle densities and, therefore, 
we neglected this effect in the theoretical model, we can- 
not completely exclude that some variations of plasma 
parameters do occur. The third source is the systematic 
difference in particle drift velocities observed when using 
different dispensers and documented for the largest par- 
ticle size investigated (see Fig. jT]). This important effect 
is under current investigation and is presumably associ- 
ated with long-living neutral flows excited by the gas-jet 
dispenser |56j . 



VII. CONCLUSION 

It has been observed that flows of particles in a positive 
column of a horizontal dc discharge are characterized by 
different velocities, depending on whether the experiment 
is performed in ground-based (laboratory) or micrograv- 
ity conditions, at otherwise identical discharge param- 
eters. In laboratory particles drift systematically faster 
than in microgravity and the velocity difference increases 
when increasing the particle size and/or decreasing neu- 
tral gas pressure. In this paper we have provided an 
explanation for this observation. Qualitatively, gravity 
shifts the particle downwards from the tube axis, where 
the radial electric field is strong enough, so that the ver- 
tical component of the electrical force can balance for the 
particle gravity. This region is characterized by faster ion 
flows and, as a result, the charge of the particles increases 
compared to the situation where the particles are located 
close to the tube axis (e.g. in microgravity). Thus, the 



particles in laboratory experiments feel stronger longi- 
tudinal electric force and hence drift faster. To be able 
to make quantitative comparison between theory and ex- 
periments we have developed an analytical model, which 
properly accounts for the main plasma-particle interac- 
tion mechanisms important for the present study. The 
proposed model yields reasonable agreement with the ex- 
perimental measurements. It can be of certain value in 
other situations characterized by relative drifts between 
the ion and particle components. 
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Appendix: Ion-neutral collision cross section in neon 

We combine the conventional definition of the effective 
ion- neutral collision frequency u ~ eE /vnii with Eq. ([5|). 
In the limit of slow drifts {u ^ wt, ) we have from Eq. ^ 
V ~ UnavTi: which results in tr ~ evxja ~ 1.3 x 10"^'* 
cm^. In the limit of fast drifts {u ^ wtJ we have v ~ 
UnCru, which results in cr ~ ev^b/a^ ~ 0.6 x 10~^^ cm^ 
(in both cases we assume Ti ^ Tn — 0.03 eV). In this 
paper we take for the cross section an "average" value 
(T ~ 1 X 10"^"* cm^, which is a good approximation for 
near-thermal drifts. 
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